Four- Wave mixing in degenerate Fermi gases: Beyond the undepleted pump 

approximation 
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We analyze the full nonlinear dynamics of the four-wave mixing between an incident beam of 
fermions and a fermionic density grating. We find that when the number of atoms in the beam 
is comparable to the number of atoms forming the grating, the dephasing of that grating, which 
normally leads to a decay of its amplitude, is suppressed. Instead, the density grating and the beam 
density exhibit large nonlinear coupled amplitude oscillations. In this case four-wave mixing can 
persist for much longer times compared to the case of negligible back-action. We also evaluate the 
efficiency of the four-wave mixing and show that it can be enhanced by producing an initial density 
grating with an amplitude that is less than the maximum value. These results indicate that efficient 
four-wave mixing in fermionic alkali gases should be experimentally observable. 

PACS numbers: 03.75.Fi,05.30.Fk 



INTRODUCTION 

In recent years the field of nonlinear atom optics 
has become firmly established both theoretically and ex- 
perimentally with the demonstration of four- wave mixing 
in Bose-Einstein condensates 0, coherent matter wave 
amplification 0], the creation of dark and bright 0] 
atomic solitons, the generation of correlated atomic pairs 
and squeezing |a,l3) and the creation of a molecular con- 
densate - the analog of second-harmonic generation in 
optics H, . So far, this work has been largely limited to 
bosonic atoms and has been based to a great extent on 
analogies with nonlinear optics: Ultracold bosons inter- 
act primarily via s-wave collisions and typically occupy 
only a few quantum states, so that their mean-field de- 
scription in terms of a Gross-Pitaevskii equation includes 
a nonlinearity of the same form as a third-order optical 
Kerr nonlinearity, and thus many of the results from non- 
linear optics can be directly applied to condensates. In 
particular, Bose-Einstein condensates are often viewed 
as a matter-wave equivalent of the optical laser- the so 
called " atom laser" . 

The last few years have also seen substantial progress 
in the cooling and trapping of atomic Fermi gases to the 
quantum degenerate regime, with temperatures reaching 
as low as T < 0.2Tp where Tp is the Fermi temper- 
ature m E [11 El- It is therefore only a mat- 
ter of time before experiments start to probe the dy- 
namics of these gases far from equilibrium and perform 
fermionic nonlinear atom optics experiments. However, 
the extension of the ideas developed in nonlinear optics 
to Fermi systems is by no means obvious or straightfor- 
ward because of the different quantum statistics obeyed 
by fermions and bosons. Whereas for bosons the anal- 
ysis can typically be limited to a few macroscopically 
occupied modes, the Pauli exclusion principle prohibits 



the number of modes from ever being less than the num- 
ber of particles for fermions. This makes both analytical 
and numerical studies of the dynamics of Fermi systems 
rather challenging, except for very limited situations (e.g. 
T = and weak perturbations). Consequently, the idea 
of doing nonlinear atom optics with fermions is still new 
and largely unexplored. 

It is only recently that it has been convincingly argued 
that four- wave mixing is possible with fermions (15» . .ld| . 
These works showed how four-wave mixing can be inter- 
preted in terms of incident particles undergoing Bragg 
scattering off of a fermionic density grating, but were 
limited to an incident beam consisting of a single test 
particle. They also neglected the dynamics of the density 
grating. Related work has shown that four-wave mixing 
between a degenerate fermion beam and a bosonic den- 
sity grating is also possible While these results are 
promising, a complete analysis of purely fermionic four- 
wave mixing is still lacking. 

In a previous paper [l^ . we studied fermionic four- 
wave mixing using the matter-wave analog of the unde- 
pleted pump approximation in optics, i.e. we neglected 
the back-action of the diffracted atoms on the density 
grating. The density grating was formed by atoms in a 
single hyperfine state prepared in a coherent superposi- 
tion of the momentum states k — q/2 and k -\- q/2 for 
all k such that |fc| < kpg, where kp^g is the Fermi wave 
number of the grating atoms and q the wave number of 
the density modulation. The third incident wave con- 
sisted of a degenerate Fermi gas in a different hyperfine 
state so that the atoms in the beam and grating could 
interact via s-wave collisions. These atoms could then 
be diffracted by the density grating to produce a fourth, 
scattered wave. Because of the range of energy states 
occupied by its constituent atoms, the amplitude of the 
grating would decay away due to the dephasing of the 



atoms. This led to a finite time during which four-wave 
mixing could occur and placed serious constraints on the 
properties of the fermions involved. 

In this paper, we extend the results of Ref. by 
taking into account the back-action of the incident beam 
on the fermionic density grating. We find that when the 
number of atoms in the incident beam becomes a siz- 
able fraction of the number of atoms forming the density 
grating, it no longer decays away due to the dephasing. 
Instead, it exhibits large nonlinear amplitude oscillations 
that are coupled to the Bragg oscillations of the beam. 
This leads to the efficient generation of the fourth scat- 
tered wave, even for times much longer than the grating 
dephasing time. We interpret this behavior in terms of 
a Bloch vector picture. We also evaluate the efficiency 
for the generation of the scattered wave and show that 
it should be possible to produce a macroscopic fourth 
wave. Interestingly enough, we find that the efficiency 
of the four- wave mixing can be improved significantly by 
initially preparing a density grating with an amplitude 
that is less than maximal. 

The organization of this paper is as follows: In section 
II, we briefly review the model developed in 0| . Section 
III presents numerical results that show how the dephas- 
ing of the grating can be eliminated when the back-action 
becomes significant. The efficiency of the four- wave mix- 
ing is also calculated and it is shown that a macroscopic 
fourth wave can be readily produced. Finally, section IV 
is a summary and conclusion. The appendix presents an 
analytically soluble model in terms of coupled Bloch vec- 
tors that yields a qualitative interpretation of our results. 

MODEL 

In this section we briefly review the one-dimcnsional 
model developed in to study four-wave mixing be- 
tween fermions at zero temperature. We consider two 
spin-polarized species of fermions- a density grating com- 
posed of spin-up polarized atoms and an incident beam 
of spin-down polarized atoms, see Fig.^ Atoms in differ- 
ent spin states interact via elastic s-wave scattering. The 
cross-section for p-wave scattering at T = is orders of 
magnitude smaller than that of s-wave scattering, thus 
p-wave scattering is neglected in our model . 

For clarity and simplicity we restrict ourselves to one 
spatial dimension. It is usually possible to reduce the 
problem of four-wave mixing to a two-dimensional pla- 
nar geometry, with one dimension parallel to the den- 
sity grating. Since the momentum transfer to the beam 
atoms is parallel to the density modulation of the grat- 
ing, the essential dynamics of the wave mixing process 
occurs only in that z-direction. We assume that the 
grating is confined between x ~ to x = d along the 
perpendicular x axis. If the average momentum, hk±, of 
the atomic beam along that direction is sufficiently large 
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FIG. 1: Model of the scattering. The spin-down polarized 
beam is scattered (s-wave) by a spin-up polarized grating. 



that k± ^ 1 /d, then the requirements of energy and mo- 
mentum conservation insure that k±_ remains practically 
unchanged by the scattering. In this case, wave mixing 
in two or three dimensions is in essence the same as in 
the one-dimensional model that we consider. (The trans- 
verse momenta of the beam atoms can at most change 
sign due to reflection off the grating at x = 17], but 
the reflected wave can be minimized by using a smooth 
grating density profile in the x-direction.) 

The second quantized Hamiltonian describing the sit- 
uation at hand is 

k 

+ (4i+9T"L-9i«fe2iafeiT) ' (1) 

with ujk = hk'^/2m and Uq ~ ATTha/{Ld^m). Here, a is 
the scattering length, d^ is the transverse cross-sectional 
area, m the mass of one atom and L the quantization 
length. The operators ak,a and a[. ^ are annihilation and 
creation operators for atoms with momentum hk and spin 
f = T? ii satisfying standard fermionic anticommutation 
relations. 

The dynamics of the system is conveniently described 
in terms of the particle-hole operators 

(b) _ t 

Pkl' = 4+fc'T"fcT. (2) 

that create a particle in state k + k' and a hole in k 
for the beam and grating, respectively. They obey the 
Heisenberg equations of motion 

^ Jb.g) _ (b,g) 
^^PkM' -^k,k'Pk,k' 

+ ^0 2^ yPk+k^,k'-kl " Pk,k'^ki) Pk2M' ^"^i 
fci ,k2 
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with u!k,k' =^k- (^k+k' ■ 

We consider a beam of Nf, atoms propagating along 
the i-direction, but with a small momentum component 
around a central value q parallel to the grating, which is 
taken to be in a superposition of the states k — q/2 and 
k + q/2 for —kp^g < k < kp^g- Here, kp^g — nNg/L is 
the Fermi momentum (in one dimension) and Ng is the 
number of atoms in the grating. Such a grating can be 
created from a stationary homogeneous gas by a two-step 
process where first a two-photon Bragg 7r-pulse imparts 
a momentum —hq/2 to each atom, followed by a second 
Bragg pulse with momentum Uq that creates the desired 
superposition (see e.g. Ref. Hq]). (Note that the band- 
width of the pulses should be greater than cok^.g-) Thus 
the initial state of the beam-grating system is 



where 



|*(0)) = ^ 

n (4-,/2t+«U/2t) |o)' 

I A; I "^k^.g 

n «uio). 

I ^ I ^^F, b 



(4) 



(5) 



The initial expectation values of the particle-hole op- 
erators p\^]., and p'"^],, are readily found to be 



p'iU0)\^^'^)=Sk'^oe{kp,t-\k-q\), 



(*^'^|pi'i'(0)|*(^)) = ^(4',o + 4'je(fc^,,-|fc+||) 

+^(4',o + 4',-,)e(fc^.,,-|fc-||). (6) 



We note that the spatial density of the Fermi gas is given 
by 



Jb,g) 



k' \ k / 



— ik'a 



(7) 



from which it is clear that the density for the initial state 
of the spin-up polarized grating atoms exhibits as ex- 
pected a spatial modulation with the grating wave num- 
ber q and an amplitude Ng/ L. In contrast, the spin-down 
polarized incident beam has a uniform density. 

It is convenient to introduce the slowly varying 
particle-hole operators p^'k) i^) = pi^'jf'^ (^)e~*'^'=''='* whose 



kp.b ~ 7riVf,/L being the Fermi momentum of the beam. Heisenberg equation of motion simply become 



k\ ,k2 



In the following we factorize the expectation values of these operators for different spin states, 



(8) 



(9) 



This ansatz results in the truncation of the BBGKY-type hierarchy of equations of motion for the higher-order 
moments of the particle- hole operators after the first moments |22j , yielding the closed set of c- number equations 



Ab.g) _ 



(10) 



A;i ,k2 



r 



where p[^l, = {p[^l) and = {p["i,). These equations 
can be readily integrated numerically. 

Implicit in Eqs. (|10(l and © are the three fundamental 
time scales that govern the dynamics of the interaction 
between the beam and the grating. The first time scale is 
associated with the kinetic energy ujq gained by an atom 
scattered by the grating. The other two time scales are 
associated with the density grating itself. If, for the mo- 
ment, we revert to the undepleted pump approximation 
used in Ref. 0|, namely p~^\{t) = pi^](0), then we ob- 



-Sb) - lf.ib) 



tain the linear equations of motion for the beam, 
d 



_ -ib) 



9{t) [Ptlq.k'-q 



Here 



-ib) 

Pk-q,k'+q 

= El 

2 

UoNg 



-ib) 

Pk,k'-q^ 



-(fc) 

Pk,k'+q^ 



(11) 



^e(fcF,g- |fc2+g/2|)e-^"'=-''* 
/ hqkpgV 



smc 



V m 
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where the second equality was obtained in the contin- 
uum Hmit, J2k ~^ L/2tt J dk. The coupUng constant g{t) 
is proportional to the amplitude of the density grating, 
which decays to zero in a time of the order of 

T]j = m/hqkp,g. 

This decay is caused by the spread in energies of the 
fcrmions that comprise the grating, HAlu = h\ujkp +q/2 ~ 
'^feFg-g/sl? so that the different energy states in the grat- 
ing get out of phase with each other over time. This leads 
to an effective grating lifetime of the order of td = 1/ Auj, 
which corresponds to the relative phases of the different 
states being of order unity. 

The third time scale is given by the characteristic time 

TB = 1/UnNg 

for atoms to scatter off of the grating in the limit of a 
static grating (td — > cx)). For an incident particle with 
momentum Tg^ is the frequency of oscillations be- 

tween the degenerate states q/2 and —q/2. In the Bragg 
regime that we consider in this paper, both the kinetic 
energy and the momentum of the atoms are conserved 
by the scattering process. This corresponds to the con- 
dition, 

WgTs > 1. 

The population that develops in states centered around 
±3g/2, ±5g/2, . . . is then at most on the order of 
(r^Wg)"^. If we choose the average momentum of the 
beam as q = q/2, with q > kpj,, only states centered 
around —q/2 are phase- matched to the states in the in- 
cident beam, so that we can restrict ourselves to those k 
states within a width 2kF,b of ±9/2. (This is in contrast 
to the Raman-Nath regime, where the change in the ki- 
netic energy of an atom as result of scattering off of the 
grating is negligible compared to the interaction energy 
with the grating HJ.) 

In the undepleted pump approximation it is important 
that the scattering time be less than the lifetime of the 
grating, 

TB <TD, (12) 

otherwise there is no significant scattering beam gener- 
ated. However, as we show in the next section, four-wave 
mixing can continue for times much longer than to when 
the effects of the back-action on the grating become sig- 
nificant. 

NUMERICAL RESULTS 

In this section we study the phenomena in fcrmionic 
four-wave mixing resulting from the back-action of the 
atomic beam on the grating. We proceed by numerically 
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FIG. 2: Populations n^''' as a function of time (in units of 
TB = 1/UoNg) for (a) Nt/Ng = 0.1 and (b) Nb/Ng = 0.5. In 
both cases td/tb = 7.5 and u)q/{UoNg)=2. 



solving Eqs. ifTUIl . subject to the initial conditions ® 
with q = q/2 > kp^t, using a fourth-order Runge-Kutta 
algorithm with fixed time steps chosen to produce a rel- 
ative error smaller than 1%. 

The back-action of the beam on the density grating is 
negligible for all times when the number of beam atoms 
is much less than the number of atoms in the grating, 
Ng ^ Ni,. The reason is that for every atom in the beam 
that suffers a momentum transfer of ihq, one atom of 
the grating sees its wave function "reduced" from the 
superposition state {\k — q/2) + \k + q/2))/\/2 to the state 
\k =p q/2). Since only those atoms in a superposition of 
k — q/2 and k + q/2 (for all possible k) contribute to 
the density modulation, the fractional reduction in the 
amplitude of the grating will be at most Nb /Ng ^ 1 and 
the undepleted pump approximation is valid for all times. 
This approximation also holds for arbitrary Nh for times 
short enough that the number of atoms scattered by the 
grating is much less than Ng, that is, for times t -^tb- 

Fcrmionic four-wave mixing is conveniently character- 
ized in terms of the occupation numbers of the momen- 
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FIG. 3: The densities \pq\ as a function of time scaled by 
UoNg for the ratios Nt/Ng = 0.1,0.5 and 1.0 in the subplot 
(a), (b) and (c), respectively. 



turn states of the beam and grating atoms 



nr'(r)-f^^'(r), (i3) 

and the amphtude of the density modulation of the grat- 
ing, 



(14) 
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FIG. 4: Atomic populations as a function of time (in units of 
tb = 1/UoNg) for the case of equal particle numbers Nb/Ng = 
1.0. The upper plots (a) show the density of the grating, n]^\ 
and the lower plots (b) show the beam evolution, ri^^^ . 



Fig. |2Jl plots rS^^ as a function of time for the case of 
a "large" grating, Nb/Ng = 0.1 and td/tb = 7.5, and 
Fig. shows the evolution of the density modulation of 
the grating |pg(t)| for these same parameters. The atoms 
in the beam undergo Bragg oscillations between the k 
states centered around —q/2 and +q/2 up until a time of 
the order of t^, when the amplitude of the oscillations 
dies out. It is clear from the behavior of |pq(i)| that the 
disappearance of the Bragg oscillations corresponds to 
the decay of the grating amplitude due to dephasing. In 
this example, the back-action of the beam on the grating 
manifests itself in the small amplitude oscillations super- 
imposed on the sine function behavior of |pg(i)|. These 
oscillations reflect the coupling between individual modes 
of the beam and of the grating, which all occur at roughly 
the same frequency t^^ for ^ |fcg/m|, as further dis- 
cussed later on. 

Figs. HJd and ^ are the same as Fig. except that 
"smaller" gratings, with Nb/Ng = 0.5 and Nb/Ng = 1.0, 
respectively, have been used. For these parameters there 
is a dramatic change in the behavior of \pq{t) \ as the back- 
action becomes more pronounced. Instead of decaying 
away, the density grating now undergoes large amplitude 



oscillations, as illustrated in Fig. O At the same time, 
the beam atoms undergo nearly perfect Bragg oscillations 
that persist for times much longer than td ■ Note that the 
amplitude of the oscillations is equal to 1 for Nb/Ng = 
0.5, but only to about 1/2 for Nb/Ng = 1.0 (see the scales 
of the grey-scale plots) although there is no evidence of 
any decay in the Bragg oscillations. 

If one ignores the weak population of the states cen- 
tered around ±3g/2,±5g/2, . . . (an approximation valid 
in the Bragg regime) these coupled oscillations involve 
the oscillation of a beam atom between the states \k' + 
q/2) and \k' — q/2), while an atom in the grating oscil- 
lates between jfc — q/2) and jfc -I- q/2), Tr/2 out of phase 
with the beam atom, as follows from the existence of the 
constants of motion 



E 



''k+q/2 ' "'k+q/2^ 

\k\<kFg |fe|<feF_i, 



E 

\k\<kF,c 



E 

\k\<kF,t, 



(15) 



(16) 



We can gain a better understanding of fermionic four- 



6 



wave mixing by recasting it in terms of the dynamics of 
an inhomogeneously broadened ensemble of coupled two- 
level systems. In this picture, an individual two-level sys- 
tem consists of the manifold of beam and grating states 
{|A: — g/2), \k + q/2)} for all states with |fc| < kp^g and 
\k\ < kp_b for the grating and beam atoms, respectively. 
We proceed by introducing the associated Bloch vectors 

El, 



Analytic Solution 
X Numerical Data 



with the components 



(9) 



n 



(a) 

k+q/2 



A9) 

''k-q/2- 



(17) 

(18) 
(19) 
(20) 



It is easily seen that the Bloch vectors associated with 
the atoms in the incident beam obey the equations of 
motion 



|ui^)(t)=ui^)(t)xRi^)(t), 



(21) 



that is, they precess around the effective field 

(i) ^UoJ2 Ui'^^ + UoY. ^k'^'y + 5^ki. (22) 



where 



6oJk = 



fikq 



(23) 



Here, J/^?-' and Vl,^ are the components of the fc'-th Bloch 
vector for the beam, U^^-*. They are defined in the same 
way as Eqs. (|17|l - (|20() for the grating. Clearly, the grat- 
ing Bloch vector U^^-* also obeys Eq. (|21|l . with the in- 
terchange h g. 

The components of the Bloch vector in the equato- 
rial {xij) plane of the Bloch sphere are a measure of the 
coherence between the momentum states \k — q/2) and 
\k -\- q/2) and therefore the sum of the projections of the 
grating Bloch vectors onto the equatorial plane is equal 
to the amplitude of the density grating. 

Initially all of the Bloch vectors for the grating point 
along the x axis, while those of the atomic beam, which 
doesn't exhibit any coherence, point toward the north 
pole (-f i) of the Bloch sphere. (That they point north 
rather than south follows from the fact that the initial 
momenta of the beam are in an interval of width kp,i, 
about q = q/2, see Eq. We have seen that the Bragg 



rib) 



regime requires additionally that q > fcF,&, so that n 

rib) 



ib) 

k-q/2 



is initially equal to zero and W;f' > 0.) 

In the absence of any coupling to the beam, Uq = 0, 
the version of Eq. H21|) for the grating shows that the 
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FIG. 5: Crosses: numerically obtained period of oscillation 
Tp of the atomic densities; Dotted line: analytic result of Eq. 
ll^ . Time in units of UoNg. 



evolution of the Bloch vectors is simply a rotation in the 
equatorial plane at frequency Suk ■ This "inhomogeneous 
broadening" causes the grating to dephase because of the 
spreading out of the Bloch vectors in the equatorial plane. 

Consider now the other extreme situation where the 
coUisional coupling between the beam and grating atoms 
completely dominates the dynamics, i.e. UqNi, ~ 

UoNg UoJ2k t^fc^V) > iJ^fe or equivalently, td tb- 
The amplitude of the effective field that drives the grat- 
ing, Uq J2k yk'^\^)i initially equal to zero since the ini- 
tial Bloch vectors of the beam point towards the north 
pole (+z) of the Bloch sphere. However, those vectors 
start rotating around the x-axis under the influence of 
the instantaneous grating effective field UoJ2k^k^\^)^ 
leading in turn to a projection along the y-axis. This 
then allows the Bloch vector of the grating to start ro- 
tating about that axis towards the north pole of the Bloch 



sphere at the instantaneous frequency Uo'^k y ^ 



vl^^ approaches the y-axis, U^^'' reaches the north pole 
and its projection on the equatorial plane vanishes and 
then changes sign. This reverses the sense of rotation of 
the beam Bloch vectors, which begin to move back to- 
ward the pole, at which point their projections on the 
y-axis likewise change sign. This reverses the direction 
of rotation of u[f , which then passes back through the 
north pole, reversing the sense of rotation of U^'''' , etc. 
We see that the beam atoms are transformed into a 
grating-like superpositions state, and at the same time 
the atoms forming the initial grating take on a state sim- 
ilar to the initial beam state. Thus, as time passes, the 
labels "beam" and "grating" become more and more in- 
terchangeable. One result of this back-and-forth is that 
the oscillation frequency of the populations is twice that 
of the coherences. Furthermore, their amplitude is only 
1/2, both for the beam and the grating, since the Bloch 
vectors are confined to the Northern hemisphere of the 



(9) 



Vj''\t). As 
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FIG. 



Atomic populations n^'''®' as a function of time (in 



units of tb ~ 1/UoNg) for the case of short dephasing time 
td/tb = 1.0, with Nb/Ng = 1.0 and iOq/NgUo = 2. The 
upper plots are for the grating and the lower plots for the 
beam. 



Bloch sphere, < wjl'''^^ < 1. This explains the reduced 
amplitude of the oscillations in Fig. 2| 

The analytic solution of this model given in the ap- 
pendix yields, for Ni, > Ng/2, an oscillation period of 



Tp = 27r/r2o = 



F((l/2 + 7Vg/4iV,)i/2,^/2), 

(24) 



where F(fc,7r/2) is the complete elliptic integral of the 
first kind 

We have numerically calculated the Fourier transform 
of nl^'^\t) for kF,b,kF,g <C q and td/tb = 2 and con- 
firmed that the populations of the grating and the beam 
oscillate at the same frequency fio independently of k. 
Fig. compares the period of oscillations Tp as a func- 
tion of iVb for fixed Ng to the analytical result of Eq. 
(|24|l . showing a very good agreement. We find that the 
dominant frequency in the evolution of Pq{t) is approxi- 
mately ^0/2, in agreement with the predictions of Bloch 
vector model. The small amplitude, high-frequency os- 
cillations in Pq(t) are due to a residual coupling to mo- 



mentum states around ±3q/2. 

The picture that we have developed indicates that 
the effects of dephasing can only be suppressed when 
td ^ Tb- For tb ^ td, the Bloch vectors of the grating 
no longer precess around the y-axis but rather around 

an axis at an angle Ok(t) = arctan(?7o X]fe' ^k^\^) I ^'^k) 
relative to the z-axis at the "generalized Rabi frequency" 

This again leads to the spread- 



. iUoT,k' "k' 

ing out of the Bloch vectors in the equatorial plane and 
hence a collapse of the density grating. This is confirmed 
in Fig. El which shows the dynamics of the populations 
forming (a) the grating and (b) the beam for N^/Ng — 1.0 
and td/tb = 1-0. It illustrates clearly the frequency de- 
pendence on k familiar from the generalized Rabi oscil- 
lations of two-level physics. 

The four-wave mixing efficiency can be quantified by 
the maximum number of beam atoms that can be scat- 
tered by the grating. Initially the beam atoms are scat- 
tered into negative momentum states, hence the peak 
efficiency is given by the first maximum of 



(25) 



\k\'ckp 



Fig. [3 shows rjmax for different dephasing times. We ob- 
serve a rapid decrease in efficiency when N^/Ng > 1/2. 
This result can be explained in terms of the Bloch vector 
model: we show in the appendix that when > Ng/2, 
the maximum number of beam atoms that can be scat- 
tered into states around ~q/2 is Ng/2. In practice, the ef- 
ficiency can exceed this maximum value due to the resid- 
ual coupling to states around ±8(7/2, which enhances 
the density modulation of the grating. Fig. [3 also il- 
lustrates the dramatic reduction in efficiency caused by 
the dephasing of the grating. Even in that case, though, 
the achievable efficiencies indicate that four- wave mixing 
should still be observable when the Fermi gases used for 
the grating and the beam contain comparable numbers 
of atoms. 

Interestingly, one can improve the efficiency signifi- 
cantly by choosing unequal probability amplitudes for 
the components of the grating, that is, by replacing the 
initial state j^f^^)) of Eq. Q by the more general super- 
position state 



\k\<kF^g 

(26) 

with c_q/2 > Cq/2 in our case. Fig. [S] plots the Bragg 
efficiency for several values of Cq/2 and C-q/2, showing a 
maximum improvement of about 25% as compared to the 
case of equal initial amplitudes, c_g/2 — C5/2 = l/\/2- 
Even though Cq/2 = c_q/2 = 1/V2 produces an initial 
density grating with the largest possible amplitude, sub- 
sequent scattering of beam atoms will reduce the size of 
the grating. In contrast, while the amplitude of the den- 
sity grating produced by an unequal initial superposition 
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FIG. 7: Diffraction efficiency r/max as a function of Ni,/Ng for 
td/tb = 7.5 (dotted- line), 2.5 (dashed line), and 1.0 (dash- 
dot line), with Uq/NgUo — 2. The solid line corresponds to 
maximum efficiency calculated from the Bloch vector model. 
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FIG. 8: Tjmax as a function of initial probability amplitudes for 
the grating |c_q/2|^ = 1- |c<,/2p, with td/tb = 1.0, Nt/Ng = 
0.8. 



state is less than optimal initially, the subsequent scat- 
tering of beam atoms from \k — q/2) to \k + q/2) increases 
Cq/2 while reducing c_g/2, resulting in a net increase over 
time in the grating amplitude. 

CONCLUSIONS 

In this paper we have investigated purely fermionic 
four-wave mixing in a manner that accounts for the back- 
action of the scattered beam on the density grating. Our 
numerical results indicate that the lifetime for fermionic 
four-wave mixing can be dramatically improved by mak- 
ing the number of atoms in the beam comparable to the 
number of atoms that form the density grating. Although 
the effects of the beam on the grating does reduce the 



overall efficiency of the process, the efficiencies that we 
have calculated indicate that it should be possible to con- 
vert in excess of 10% of the atoms in the incident beam 
into a fourth scattered beam. 

It is well known that four-wave mixing with bosons 
can lead to the generation of squeezed states [3 US- 
Therefore a central feature of future work will examine 
the statistical properties of the scattered beam generated 
by the four- wave mixing to see if it is possible to generate 
novel fermionic states. This will involve going beyond the 
factorization ansatz in Eq. Q to calculate higher order 
moments for the particle hole operators. Phase conju- 
gation should also be possible if one uses a superfluid 
Fermi gas due to the presence of anomalous moments, 
(afc|a_fc|), resulting from the formation of Cooper pairs. 
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APPENDIX 

In this appendix we solve the coupled Bloch equations 
for the beam and grating with dcuk = 0. By introducing 
the normalized vectors 

k 

y(b.9) ^J2vl:''''^/N,,g, (28) 
fc 

yji<^'9) =J2wt''''^/m,g, (29) 
k 

we obtain the equations of motion 

v(^) = UoNgW^'^^u^s)^ (31) 

= J7o^3(u(^)w(f) - w('')m(9)), (32) 

li^s) = -UoNbW^s)y{b)^ (33) 

{,(9) = UQNbW^s)u{b)^ (34) 

= L/oA"b('"'-^-'w'-^-' — w'-^-'w^''''), (35) 

subject to the initial conditions ^^^^(0) — w^''\0) = 1. 
Here, the ' represents differentiation with respect to time. 
We note that the constants of motion given by (|15|l and 
()16|l can be expressed as 

Nbw'-''^ + iVgW^f) ^ Nb. (36) 

For our particular initial conditions we have u'-''\t) = 
v(9)(t) = for all times. This allows us to parameterize 
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the solution in terms of the angles 9b{t) and Og{t), of the 
beam and grating Bloch vectors relative to the z-axis in 
the yz- and xz-planes, respectively, 



cos 9b{t), 
sin 0b{t), 

cos Og{t), 

sin Og{t). 



(37) 
(38) 
(39) 
(40) 



By introducing the sum and difference angles 

e± = 0g ± 06 - TT, 

we obtain the uncoupled equations 

e± + u? sin e± = 0, 



(41) 



where u: — Uq ^NgN^- Eg. (jlTT) is the equation of motion 
for a plane pendulum 26] . 

Since the initial conditions for these pendula are 
^±(0) = -vr/2 and ^±(0) = ±NgUQ, a small amplitude 
linearization of the equations of motion is not justified. 
However it is still possible to find a solution. We can 
express the total "energy" for the pendula as 



sin2(6'±/2) = 2u?K^ 



where 



1 

4 7Vf,' 



and d±fi — 2arcsinK arc the turning points for the pen- 
dula. If A'';, < Ng/2, then there are no turning points and 
the pendula undergo complete revolutions about "their" 
axes. In this case the Bloch vector for the beam under- 
goes complete revolutions in the yz-plane as can be seen 
fromEq. PHji . 



If, on the other hand, 
oscillatory solution 27], 



Nfj > Ng/2, then we have the 



sm 



-Xt)/2) = ±Ksn{uj{tTto),K) 



(42) 



where sn is a jacobian elliptic function 28] and to is de- 
termined by the initial condition, Ksn(ujto , k) — 1/V2. 
In this case the beam Bloch vector has a turning point 
at cos 9b = 1 — Ng/Nb, so that the maximum number 
of beam atoms that can be scattered into states around 
— g/2 is Ng/2. The period of the oscillations for the el- 
liptic function is given by Eq. 1)24(1 . 
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